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ABSTRACT 

A method is proposed for optimizing the fabrication process of unidirectional metal- 
matrix composites. The temperature and pressure histories are optimized such that the 
residual microstresses of the composite at the end of the fabrication process are minimized 
and the material integrity throughout the process is ensured. The response of the com- 
posite during the fabrication is simulated based on a nonlinear micromechanics theory. 
The optimal fabrication problem is formulated and solved with non-linear programming. 
Application cases regarding the optimization of the fabrication cool-down phases of unidi- 
rectional ultra-high modulus graphite/copper and silicon carbide/titanium composites are 
presented. 
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Young’s Modulus. 
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Pressure. 

Strength. 
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1. INTRODUCTION 


During the last few years, significant interest has developed for metal-matrix com- 
posite (MMC) materials. Despite their high cost, metal- matrix composites (MMCs) are 
promising candidate materials for applications demanding high operational temperatures 
(from 800°F to 2000°F), high specific stiffness and strength, hygrothermal resistance, sta- 
bility, and tailorable mechanical properties. A problem with most MMCs is the high 
residual thermal microstresses developed during the cooling phase of the fabrication pro- 
cess. These high residual stresses are primarily caused by the large temperature differential 
imposed during the cooling phase because of the mismatch between the thermal expansion 
coefficients of the composite’s constituents. 

The development of thermal residual stresses during fabrication is a known problem 
for composite materials. Much research has been reported on the thermal residual mi- 
crostresses of polymer-matrix composites. We selectively mention some early studies on 
thermal stresses of composite laminates [1-2], as well as, some studies on the prediction of 
residual stresses during the curing of thermoset-matrix composites [3] and during the solid- 
ification of thermoplastic-matrix composites [4-6]. However, the residual thermal stresses 
generally are more critical for MMCs than polymer matrix composites because of the much 
larger temperature differential involved in the fabrication of MMC’s. 

The residual microstresses degrade the mechanical properties of MMCs, may initiate 
matrix failures, and are primarily responsible for their poor thermo-mechanical fatigue 
endurance. Experimental and computational studies [7-8] indicate that in addition to 
other factors, the residual microstresses depend on the actual cool-down and consolidation 
histories during fabrication, since the metallic matrix undergoes in situ nonlinear deforma- 
tions. The material behavior is further complicated by the fact that the in situ mechanical 
properties of the matrix depend on temperature and stress. Hence, a process may exist 
with optimal temperature and pressure variation, that will produce MMCs with favorable 
thermal residual stresses and improved thermo-mechanical fatigue (TMF) life. 

To the authors’ best knowledge, no research has been reported on the prediction of 
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optimal cooling and consolidation histories for MMCs. Research on the optimal fabrication 
process of polymer-matrix composites [9] has limited applicability because in most cases 
the highly non-linear matrix behavior is not an issue. The present paper describes recent 
studies on the development of a methodology for optimizing the fabrication parameters of 
unidirectional MMCs such that: (1) critical residual microstresses are minimized; and (2) 
failures in the matrix, fibers, and fiber/matrix interphase are prevented during fabrication. 
All these are directly related to TMF endurance and ensure its improvement, while item 
(2) also ensures the integrity of the composite. 

The residual microstresses and the mechanical properties of the material are predicted 
based on integrated non-linear micromechanics for metal-matrix composites developed at 
the Lewis Research Center [10-11]. The optimization problem is formulated and solved 
with the modified feasible directions mathematical programming method. The material 
nonlinearity requires special considerations in the formulation and solution of the problem 
which are addressed herein. Applications of the developed methodology are performed on 
the optimization of the cool-down process of unidirectional silicon-carbide(SiC)/titanium 
(Til5-3-3-3) and ultra high modulus graphite(P100)/copper composites. 

2. FABRICATION PROCESS 

The variety in properties of fibers and matrices, the size of the fibers, the chemical 
reactivity between fibers and matrix at higher temperatures, and the attained quality and 
properties of the final composite have resulted in the development of various fabrication 
techniques. Among them, superplastic diffusion bonding and hot isostatic pressing appear 
to have many merits and potential. A typical fabrication cycle, as the one shown in Fig. 1 
for P100/ copper, usually involves three phases. During Phase 1, the temperature of the raw 
materials is elevated near the matrix melting temperature. Application of consolidation 
pressure is usually required in Phase 2 such that the matrix is diffused and bonded with 
the fibers. It is assumed that Phase 2 is sufficiently long such that perfect bonding is 
accomplished and consolidation stresses are negligible. The final cool-down phase (Phase 
3) follows in order to reduce the temperature and pressure to the reference conditions. 
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Phase 3 seems to be most critical in the subsequent performance of the fabricated 
MMC, in that the residual microstresses are substantially developed during this phase as 
a result of the mismatch between the fiber and matrix thermal expansion coefficients. Me- 
chanical microstresses are also present as a result of the consolidation pressure. Part of the 
mechanical microstresses vanishes when the consolidation pressure is removed, however, 
their simultaneous action with the thermally induced microstresses may contribute to ma- 
trix failures. Since the development of residual stresses and the integrity of the composite 
material is primarily affected by the temperature and pressure histories in Phase 3, it is 
reasonable to focus the current study on the cool-down phase. 


3. THERMO-MECHANICAL RESPONSE OF MMCs 
The simulation of the thermo-mechanical response of MMCs during Phase 3 of the 
fabrication process requires nonlinear micromechanics capable of handling, among other 
factors, the effects of temperature, the nonlinearity of the constituent materials stress- 
strain behavior and the residual microstress build-up. Such an integrated micromechanics 
theory developed by Chamis and co-workers [10-11] is utilized in the present study. In 
summary, the theory proposes the following relation of in situ constituent properties P to 
state- variables such as temperature, stress, and time: 


p_ _ r rw ~ T S - 0- <F - t , p 

Po [ T m - r o J [ So i [ t F ] 


(1) 


The exponents in the previous equation may be estimated from experimental data. Can- 
didate properties for this equation are the moduli, Poisson’s ratios, strengths, and thermal 
expansion coefficients. 

The micromechanics theory assumes the composite microstructure shown in Fig. 2 
which consists of three material regions: the fiber, the matrix, and an interphase between 
them. The interphase represents either a special fiber coating or other substances such as 
products of chemical reactions between the fiber and matrix. Three distinct microregions 
may be recognized in the composite material. The three regions are shown in Fig. 2 and 
are identified with letters A (matrix), B (matrix-interphase), and C (matrix-interphase- 
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fibers). Externally applied mechanical and thermal loads on the composite will result in 
the development of different average stresses in the three microregions [10-11]. Because 
of the material non-linearity expressed by eq. 1, the calculation of composite properties 
and microstresses at each time step of the simulated fabrication phase requires an iterative 
solution of the governing equations. 

4. OPTIMAL FABRICATION PROCESS 

As stated previously, it is desireable to optimize the process parameters (temperature, 
pressure, and time) during the cool-down phase of the fabrication cycle such that: (1) the 
residual microstresses in the fabricated composite are favorably controlled, and (2) the 
integrity of the resulting composite is ensured. Considering the large number of parameters 
and the complexity of the simulation, this may be best accomplished with non-linear 
programming (NLP). A standard constrained NLP problem has the form: 

minimize F(z) 

subjected to constraints: 


L < z < z u 

(2.2) 

Q( z) < 0 

(2.3) 

H( z) = 0 

(2.4) 


where F( z) is the objective function, Q{ z) is the inequality constraint, H{ z) is the equality 
constraint, and z refers to a set of design variables. For our case, the design variables consist 
of a set of temperatures, pressures, and times sufficiently defining the cool-down phase. 

The desirable state of residual microstresses is related to the projected service require- 
ments of the composite. In the rest of the paper, the authors assume that reduction of 
residual microstresses is most desirable. In this context, the optimal fabrication problem 
is first formulated as the following constrained minimization: 

min(max{u?i<T^ 11 , W2&m22 ) W 3 a 'm22 > w 4 a rn22}) ( 3 ) 


( 2 . 1 ) 
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subjected to fiber and matrix stress failure constraints at n s time steps throughout the 
solution. 


iSfnch < ( a fn)j < (SfuT)j 

j = l,...,n # 

(4.1) 

{Sf22c)j < ( cr /22); < f22T )j 

£ 

f-H 

II 

' ^ 

(4.2) 

(S£nc)i < (4i); < 

j = l,...,n s 

(5.1) 

(^m22c); < i a m22)j < i^m22T)j 

j = 1, 

(5.2) 

{Sm22c)j < ( a m22)j < ($ m22r)j 

II 

1— ' 

(5.3) 

{S m 22 c)j < (°’m 22 ); < {^m. 22 T)j 

j = l,...,n. 

(5.4) 


Upper and lower bounds on the optimization vector z are also imposed in accordance with 
eq. 2.2. 

Although eqs. 3-5 adequately represent the stated objectives, a difficulty immedi- 
ately arises when a microstress increases beyond the ultimate strength value. In such a 
case, the material fails resulting in elimination of the load caring capacity of the respec- 
tive microregion. Therefore, the previously formulated constrained optimization problem 
is only meaningful within the feasible domain. Moreover, the objective function and the 
constraints are continuous in the feasible domain only since discontinuities occur during 
each failure. The discontinuity does not preclude the existence of a local minimum in the 
feasible domain, however, the matrix nonlinearity (eq. 1) mandates that non-linear pro- 
graming algorithms designed to operate outside the feasible domain (eg. exterior penalty 
methods) are inapplicable to the present problem. In addition, the starting optimization 
process should be feasible, i.e. no microfailures should occur during fabrication. The last 
requirement would eliminate many of the merits of the proposed formulation in eqs. 3-5, 
since the prediction of an initial fabrication process without material failures may be cum- 
bersome. In order to overcome these difficulties, the final temperature is included in the 
objective function (eq. 3): 


T 

rnin(max{w 1 cr^ 11 ,W2^ 22 ,W3^m22 5 u; 4 0-m22} + w 57f-) 


( 6 ) 
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The inclusion of final temperature in the objective function allows the use of a nonambient 
final temperature in the starting process. The prediction of an initial fabrication process 
with a sufficiently high final temperature such that microfailures do not occur is trivial. If a 
feasible fabrication process exists with an ambient final temperature, then the minimization 
of the proposed objective function will reduce the final temperature to the room value. In 
this manner, the associated difficulties with the knowledge of an initial feasible process are 
overcome and the use of the current methodology may be extended in the prediction of 
feasible fabrication processes. 

The optimization problem described by eqs. 6 and 3-5 is numerically solved with 
the modified feasible directions NLP method [12]. One important feature of the modified 
feasible directions algorithm is that it performs a direct search within the feasible opti- 
mization domain based on the gradients of the objective function and the constraints. The 
algorithm includes also an active set strategy, i.e. only the constraints near violation are 
included in the estimation of the search direction enabling the efficient handling of the 
large number of constraints defined in eqs. 4-5. 

5. APPLICATIONS AND DISCUSSION 

This method was applied to optimize the cool-down phase of the following two 
unidirectional MMCs: (1) 0.40 fiber volume ratio (FVR) ultra high-modulus graphite 
(P100)/copper, and (2) 0.40 FVR SiC/Til5-3-3-3. Fabrication data regarding the current 
fabrication processes of these composites were provided by the Materials Division of NASA 
Lewis Research Center. Representative constituent properties of both composite systems 
at reference conditions (70°F, 0 psi) are shown in Tables 1 and 2. 

For the analysis, phase 3 was subdived into four increments of linearly varying tem- 
perature and pressure. Stress constraints were imposed at five evenly spaced time inter- 
vals in each linear segment. In this manner, twenty constraints were introduced for each 
microstress inequality described in eqs. 4-5. Only the residual microstresses and 

<r ^22 were included in the objective function (eq. 6), that is: w i = u >2 ~ w 5 = 1 and 
U 7 3 = — 0. The temperatures, pressures, and times at the starting and final points of 
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the four linear segments were used as optimization parameters. The temperature at the 
beginning of Phase 3 was held constant and equal to the respective temperature of the 
current processes, and the final pressure was set equal to zero. The following upper and 
lower bounds were imposed on the optimization variables in accordance with eq. 2.2: 


T 0 <T<T m 

(7.1) 

0 < p < 50 ksi 

(7.2) 

lOsec < t < 18000.sec 

(7.3) 


Case 1: PlOO/Copper 

Fig. 3 shows the current and the resultant optimum fabrication process for the 
PlOO/copper MMC. The cool-down phase of the current process consists of two segments 
of linearly decreasing temperature and pressure. In the resultant optimum process the 
temperature remains high during the first 20Q0sec, is reduced to room temperature in the 
following 2500sec, and remains there until the end of the process. The predicted optimal 
consolidation pressure in Fig. 3b is significantly higher than the pressure of the current pro- 
cess. At Osec the pressure is 1.58ksi, in the next 4500sec the pressure gradually increases 
to 6.5ksi, and then drops to zero. The pressure increases as the temperature decreases 
makes sense, because in this manner the matrix becomes highly nonlinear and undergoes 
“plastic” deformation resulting in reduced matrix microstress ll . The predicted build- 
ups of microstresses <7^ u and 2 are shown in Fig. 4. The final microstresses in Fig. 4 
are the residual microstresses in the fabricated composite. The residual microstress cr^ n 
of the optimal process decreased by 14.5% compared to the respective microstress value of 
the current process. The residual microstress cr^ 22 of the optimal process was increased 
and both residual microstresses of the optimum process have nearly equal values. The 
predicted microstress build-up in the current process is more uniform compared to the op- 
timum process because both pressure and temperature decrease almost linearly, as opposed 
to the more complicated temperature-pressure variation of the optimal process. Figs. 5 
and 6 illustrate the variation of the longitudinal and transverse in situ matrix moduli and 
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strengths respectively during the current and the optimal process. The fact that the ma- 
trix in the optimal process is undergoing nonlinear deformations is further illustrated in 
Fig. 5 by the variation of the in situ moduli. The high residual microstresses resulted 
in a significant decrease of the matrix longitudinal modulus in both current and optimal 
processes. The optimal process has a higher final in situ matrix longitudinal modulus, and 
decreased final transverse modulus compared to the current process. The in situ strengths 
of the matrix did not change. 

Case 2: SiC/Til5-3-3-3 

The resultant optimum process for SiC/Til5-3-3-3 is shown in Fig. 7 together with 
the current one. The microstress build-up is shown in Fig. 8, and the variations of the 
in situ moduli and strengths are respectively shown in Figs. 9 and 10. In contrast to the 
previous case, the process optimization did not result in any significant reductions of the 
final microstresses with respect to the currently used process. A possible explanation is 
that the current process may have been improved through a trial-and-error procedure. 

As previously stated, the principal physical mechanism that has contributed to mi- 
crostress reductions in the previous case studies was the non-linear in situ matrix behavior, 
as a result of high temperature and pressure. Two additional physical phenomena not in- 
cluded in the present study, namely the matrix solidification history and viscoplasticity, af- 
fect the microstress development in the composite. During the first stages of the cool-down 
phase when the temperature and pressure are high, the matrix actually remains in a near 
fluid state, and temperature and pressure variations do not contribute to the development 
of residual microstresses. In fact, the resultant optimal process for the graphite/copper 
composite indicates that the matrix moduli should remain low during the first stages of the 
process and illustrates the need for incorporating the matrix solidification history during 
the fabrication. Prolonged matrix nonlinearity will also result in partial relaxation of the 
induced thermal matrix microstresses. The effects of matrix solidification history require 
further consideration. Work is currently under development to include this effect in the 
simulation and optimization of the fabrication process. 
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6. SUMMARY 


A method is proposed for optimizing the fabrication process of unidirectional metal- 
matrix composites. The response of the fabricated MMC was simulated based on nonlinear 
micromechanics. The temperature and pressure variation during the cool-dowm phase of 
the fabrication process was optimized for reduced residual microstresses in the fabricated 
composite and elimination of failures. An in-house research code has been developed 
incorporating this method. 

Case studies were performed on ultra-high modulus graphite (P100)/copper, and sil- 
icon carbide (SiC)/titanium (T515-3-3-3) composites. The predicted optimal process for 
graphite/copper resulted in an estimated reduction of the maximum final microstress by 
15%, by favorably optimizing the nonlinear in situ matrix behavior. The optimization of 
the fabrication process for SiC/Til5-3-3-3 produced insignificant reductions in the residual 
microstresses. 
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Table 1 . Representative constituent mechanical properties of P100 Graphite/Copper 

at reference conditions. 


P100 Graphite Copper 


Efn — 105.0 Mpsi E m — 17.7 Mpsi 

E/22 = 0.90 Mpsi 

G fi2 = 1-10 Mpsi G m = 6.81 Mpsi 

(7/23 = 0.70 Mpsi 

pf = 0.078 Ib/in 3 p m = 0.32 Ibjin 3 

Ufi 2 = 0.20 injin = 0.30 injin 

uj 23 = 0.25 injin 

Q /n = -0.90 pin/in j°F a m = 9.80 pinjinj 0 F 

a /22 = 5.60 pinjinj°F 

Sfn,T = 325.0 ksi S m n — 32.0 ksi 

Sjn t c — 200.0 ksi 
5/22 = 25.0 ksi 

5/12 = 25.0 ksi S ms = 19.0 ksi 

5/23 = 12.5 ksi 


Table 2. Representative constituent mechanical properties of SiC /Til 5-3-3-3 

at reference conditions. 



E/n = 62.0 Mpsi E m = 12.3 Mpsi 

E j 22 = 62.0 Mpsi 

G 712 = 23.8 Mpsi Gm = 4.659 Mpsi 

<7/23 = 23.8 Mpsi 

pj = 0.11 lb jin 3 Pm = 0.172 Ibjin 3 

u fi2 = 0.30 in jin i/ m = 0.32 in jin 

v f 2 z — 0.30 injin 

a fn = 1-80 pinjin/°F a m = 4.50 pinjinj°F 

a f 22 — 1-80 pinjinj°F 

5/n, r = 500.0 ksi S mn = 130.0 ksi 

5/n, c = 650.0 ksi 

5 / 22 , T = 500.0 ksi 

5/22,c = 650.0 ksi 

5/12 = 300.0 ksi Sms = 90.0 ksi 

5/23 — 300.0 ksi 
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(a) TEMPERATURE. (b) CONSOLIDATION PRESSURE. 

FIGURE 3. - OPTIMUM AND CURRENT COOL-DOWN PHASES FOR 0.40 FVR P100/C0PPER (CASF 1). 



TIME, sec 

riGIJRF 'l. MATRIX MICROSTRESSES DEVELOPED DURING THE 
COOLDOWN PHASE OF 0.40 FVR P100/C0PPER (CASE 1). 
OPTIMUM AND CURRENT PROCESSES. 



TIME, sec 


FIGURE 5. - IN SITU MATRIX MODULI DURING THE COOL-DOWN 
PHASE OF 0.40 FVR P100/C0PPER (CASE 1). OPTIMUM AND 
CURRENT PROCESSES. 
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TIME, SEC 

FIGURE 6. - IN SITU MATRIX STRENGTHS DURING THE COOL- 
DOWN PHASE OF 0.40 FVR P100/C0PPER (CASE 1). OPTIMUM 
AND CURRENT PROCESSES. 



(a) TEMPERATURE. (b) CONSOLIDATION PRESSURE. 

FIGURE 7. - OPTIMUM AND CURRENT COOL-DOWN PHASES FOR 0.40 FVR SiC/T i 15-3-3-3 (CASE 2). 




KSI 



TIME, sec 

FIGURE 8. - MATRIX MICROSTRESSES DEVELOPED DURING THE 
COOL-DOWN PHASE OF 0.40 FVR SiC/Ti 15-33-3 (CASE 2). 
OPTIMUM AND CURRENT PROCESSES. 



FIGURE 9, - IN SITU MATRIX MODULI DURING THE COOl -DOWN 
PHASE OF 0.40 FVR SiC/Ti 15-3-3-3 (CASE 2). OPTIMUM AND 
CURRENT PROCESSES. 



FIGURE 10. - IN SITU MATRIX STRENGTHS DURING THE COOL- 
DOWN PHASE OF 0.40 FVR SiC/Ti 15-3-3-3 (CASE 2). OPTI- 
MUM AND CURRENT PROCESSES. 


17 


!M?SLnd Report Documentation Page 

Space Administration 

1. Report No. 2. Government Accession No. 

NASA TM- 102559 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

Optimal Fabrication Processes for Unidirectional Metal-Matrix Composites: 
A Computational Simulation 

5. Report Date 

6. Performing Organization Code 

7. Author(s) 

D.A. Saravanos, P.L.N. Murthy, and M. Morel 

8. Performing Organization Report No. 

E-5379 

10. Work Unit No. 
5 10-0 1-0 A 

9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Prepared for the 35th International SAMPE Symposium and Exhibition, Anaheim, California, April 2-5, 1990, 
D.A. Saravanos, Case Western Reserve University, Cleveland, Ohio 44106 and NASA Resident Research 
Associate at Lewis Research Center. P.L.N. Murthy, NASA Lewis Research Center. M. Morel, Sverdrup 
Technology, Inc., NASA Lewis Research Center Group, Cleveland, Ohio 44135. 


V 


16. Abstract 


^ . 


A method is proposed for optimizing the fabrication process of unidirectional metal-matrix composites. The 
temperature and pressure histories are optimized such that the residual microstresses of the composite at the end 
of the fabrication process are minimized and the material integrity throughout the process is ensured. The 
response of the composite during the fabrication is simulated based on a nonlinear micromechanics theory. The 
optimal fabrication problem is formulated and solved with non-linear programming. Application cases regarding 
the optimization of the fabrication cool-down phases of unidirectional ultra-high modulus graphite/copper and 
silicon carbide/titanium composites are presented. 


17. Key Words (Suggested by Author(s)) 

Metal-matrix composites; Fabrication; Optimization; 
Optimal temperature; Optimal consolidation pressure; 
Residual stresses; Nonlinear micromechanics 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 24 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of pages 

Unclassified 

Unclassified 

18 


22. Price* 

A03 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 




